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ABSTRACT 

The IfA Deep survey uncovered ~130 thermonuclear supernovae (TNSNe, i.e. Type la) 
candidates at rcdshifts from z = 0.1 out to beyond 2=1. The TNSN explosion rates derived 
from these data have been controversial, conflicting with evidence emerging from other surveys. 
This work revisits the IfA Deep survey to re-evaluate the photometric evidence. Applying the 
SOFT program to the light curves of all SN candidates, we derive new classification grades and 
rcdshift estimates. We find a volumetric rate for z ~ 0.5 that is substantially smaller than the 
originally published values, bringing the revised IfA Deep rate into good agreement with other 
surveys. With our improved photometric analysis techniques, we are able to confidently extend 
the rate measurements to higher rcdshifts, and wc find a steadily increasing TN SN rate, with no 
indication of a peak out to z = 1.05. 



1. Introduction 

In 2001, the IfA Deep Survey was undertaken as 
a collaborative program of astronomers at the Uni- 
versity of Hawaii's Institute for Astronomy. The 
survey was designed with a strong emphasis on 
discovering high redshift (z > 0.5) SNe. Over 
a 5 month span the survey produced 23 spectro- 
scopically confirmed high-z thermonuclear super- 
novae (TNSNe, or Type la), doubling the sample 
at z > 0.7 and providing an import ant verifica- 



tion of the evidence for dark energy (IBarris et al 
20041 ). In addition, a larger sample of 133 pho- 



tometrically classified SNe were used to measure 



the h igh-z TN SN explosion rate (jBarris fc Tonrv 



20061 hereafter BT06). Subsequent measurements 
from more recent surveys have consistently found 
lower values for the TN SN rate a t z ~ 0.5 than the 
BT06I results dNeill et al.l 120061 : iBotticella et~al 



20081 : iKuznetsova et all l2008h . The 



common ex- 
planation for this disagreement has been that the 



BT06I rate is affected by systematic errors in the 
photometric classification and rcdshift estimation 
- which are not present in the spectroscopic sam- 
ples of later surveys. 

Resolving this discrepancy is an important step 
toward preparing for the next generation of very 
large scale SN surveys. Spectroscopic classifica- 
tion and redshift determination are becoming un- 
realistic for all-sky survey programs yielding thou- 
sands of SNe per year. In order to make full use of 
these new kilo-SN data sets we must have a high 
confidence in our photometric analysis methods so 
that truly divergent results are not mistakenly as- 
cribed to systematic biases. 

This work describes a comprehensive reanaly- 
sis of the IfA Deep Survey data using Supernove 
Ontology with Fuzzy Templates (SOFT), a fuzzy- 
set-based method for SN classification and redshift 
estimation th at does not require any s pectr oscopic 
information (jRodnev fc Tonrvi 120091 . boioh . The 



1 



end result of this analysis is a recalculation of the 
TN SN rate as a function of redshift, extending it 
out to z=l.l. 



ment, and argued that the most likely source of 
the discrepancy is low-redshift CC SNe that were 
misidcntified as TN SN. 



The 2006 Rate Estimate 



3. SOFT Reanalysis 



The IfA Deep Survey was the first survey to 
take full advantage of the wide field imaging cam- 
eras of Mauna Kea to allow for both SN discovery 
and simultaneous follow-up observations. Using 
the 12k camera on the Canada France Hawaii Tele- 
scope (CFHT) and Suprime-Cam on Subaru, the 
survey covered 5 fields with a total survey area of 
2.5 deg 2 . Each field was revisited approximately 
every 10 days for observations in the R, I and Z 
bands. The survey design and operations ar e de- 
scribed in greate r detail bv lBarris et al.l (|2004f ) and 



Liu et al 



(120021) . 

The most promising SN candidates were quickly 
observed spectroscopically to confirm them as SNe 
and to measure their redshifts. This produced the 
spectroscopi c samp le of 23 TNSNe presented by 



Barris et al.l ( 20041) . After the conclusion of the 



5 month survey, the images were reprocessed to 
include all available photometric d ata in a sin- 
gle a nalysis. The NN2 algorithm (jBarris et al 



20051 ) was used to construct light curves from dif- 



ference i mages for over 10,000 possible transient 
objects. IBT06I applied a succession of screening 
criteria to identify 133 objects as possible TN SN 
candidates. Q This list was reduced to 98 using 
the original Baycsian Adaptive Template Match- 
ing code (a precursor to the SOFT program used 
in this work). It was this final collection of 98 
photometrically classified TNSNe that was used 
to determine the IfA Deep SN rates. 

BT06I recognized that the IfA Deep rate esti- 



mates were not in line with other published mea- 
surements. In particular, the rate at z=0.55 is dis- 
crepant by mo re than 3<7 when compared to the 

~ (l2002h. Subsequently, 



z=0.55 rate of iPain et al 



the principal critique of the BT06I rate estimate 
has been the possibility of a systematic bias due 
to misclassification of core collapse ( CC) SNe as 
TN SN . This point was emphasized by iNeill et al.l 
(|2006l ). presenting the TNSN rate measurements 
of the SNLS project. They found a rate at z=0.5 
that is > 5cr less than the lBT06l rate measure- 



1 Note that in iBarrid l|2004T l and iBTOfj there are actually 
only 131 unique objects, as two SNe were counted twice 
due to a mistake in tabulation. 



The evolution of the light curve analysis code 
from BATM in BT06 to SOFT in this work has 
introduced three significant improvements: 

• The BATM program used a set of 20 TN SN 
light curve templates with no CC templates, 
SOFT uses 26 TNSNe and 16 CC SNe. 

• BATM uses a purely Bayesian formalism, 
while SOFT uses fuzzy set theory to account 
for the incompleteness of the set of available 
templates. 

• The BATM code as applied in BT06 pro- 
vided a binary classification (la or not) and 
assigned a single redshift to each object. 
SOFT gives a classification probability and 
a redshift probability distribution. 

These advances give SOFT much greater lever- 
age for distingu i shing TNSNe from CCSNe 



(|Rodnev fc Tonrvl 120091) and for estimating the 



redsh ift from the light curve alone (jRodnev &: Tonrv 
20101) . 



Applying SOFT to each o f the 1 3 1 SN candi- 
date light curves published in lBarris ( 2004 ) yields 
a classification grade 7tn, which approximately 
reflects the likelihood that each candidate belongs 
in that classH For priors on the location parame- 
ters (time of peak t pk , host galaxy extinction Ay, 
redshift z, and distance modulus /i e ), we used the 
same "uninformative priors" that were applied in 



the validation sets in Rodney fc Tonrv ( 2009t ) 



except when a redshift measurement for the SN or 
its host is available. 



For the 23 TNSNe from iBarris et al.l (|2004f ). 



which have precise spectroscopic redshift measure- 
ments, we applied the SOFT analysis twice. The 
first pass used only the weak, uninformative prior, 
while the second pass utilized a much stronger 
spectroscopic redshift prior. This subset provides 
us with a useful internal consistency check, to val- 
idate the SOFT redshift estimates of the remain- 
ing 98 SNe. Table Q] reports the results of the 
first pass, ignoring the spectroscopic information. 



Also called the Posterior Membership Grade, PMGq 
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Fig. 1. — Comparison of the peak redshift esti- 
mates from SOFT against the spectroscopic red- 
shift measurements for 23 TN SNe. The straight 
line shows perfect parity (zsOFT=Zspec)- The rms 
scatter about the line is 0.10, and the reduced x 2 
statistic is 0.85. 

The classification grade 7tn is reported for each of 
the 23 spectroscopically confirmed TN SNe, along 
with the measured spectroscopic redshift and the 
SOFT redshift estimate when SOFT uses the un- 
informativc z prior. 

It is encouraging to note that all 23 spectro- 
scopically confirmed TNSNe are correctly classi- 
fied by SOFT (the lowest 7tn classification grade 
is around 64%). This suggests a low incidence of 
false negatives (mislabeling a TN SN as a CC SN), 
although this should not be taken too far, as 
one would expect that the sub-sample selected 
for spectroscopic follow-up includes the most eas- 
ily identified TNSN candidates. Unfortunately 
the spectroscopic sample cannot tell us anything 
about the rate of false positives (classifying a 
CCSN as a TNSN). The comparison of z SO ft 
against z spoc is shown in Figure [TJ The mean 
zsoft — ^spec residual is 0.024, and the rms scatter 
of the residuals is 0.10, with a reduced \ 2 statistic 
of 0.85. 

The SOFT classification grades and redshift es- 
timates for all 131 TNSN candidates are summa- 
rized in Table [5J For each SN in this table we 



have used the best available redshift measurement 
to inform the prior. The 23 spectroscopically con- 
firmed SNe from Table Q] appear again in Table [2j 
now with updated classification grades and im- 
proved parameter estimates due to the addition of 
spectroscopic information. Also included are red- 
shift priors derived from 14 likely SN host galaxies 
and photometric redshift estimates for a further 
22 possible hosts (see the final two columns of Ta- 
ble [2]). To account for possibile host galaxy mis- 
matches, the floor of the Gaussian redshift prior 
has been raised to 10% of the peak value, spread 
across the width of the redshift parameter space. 

The majority of these 131 candidates are clas- 
sified by SOFT as TNSNe: 90 objects have a 
TN SN classification grade above 50%, with 63 of 
those above 90%. This is, however, somewhat 
lower than t he 98 objects identified as TNSNe 
m the IBT06I analysis. As will be explained be- 
low, we include all 131 objects in our analysis, 
using the classification grades to account for un- 
certainty in sample selection. The mean redshift 
uncertainty determined by SOFT for the complete 
set is 0.055, broadly consistent with the scatter ob- 
served in FigurcQ]and seen in the SOFT validation 
tests using SNe from the Sloan Digital Sky Survey 
(SPSS) and the Superno va Legacy Survey (SNLS) 
(|Rodnev fc Tonry||2010l) . 

A handful of SN candidates in Table[2]show dis- 
crepancies between the redshift prior and the red- 
shift estimate determined by SOFT. Three objects 
(f0848d, f0848f, f0848h) have marginally signifi- 
cant (< 2a) differences between the assumed host 
galaxy redshift and the SOFT redshift. These host 
galaxy redshift priors are ta ken from the SPSS 



photometric redshift catalog (jOvaizu et al.l 120081 ) 



and have very large uncertainties (er 2 > 0.1), so 
we accept the SOFT distribution as a reasonable 
estimate. Two other objects (SN2001ft, f0848b) 
have SOFT redshifts that are significantly lower 
than the input redshift priors. Both of these 
objects have a classification grade of 7tn=0, in- 
dicating that SOFT strongly classifies them as 
CC SNe. In such cases the SOFT redshift estima- 
tion is very poor, because CC SNe are an inhomo- 
geneous population that is only sparsely sampled 
by the SOFT template library. Our 7TN-weighted 
counting procedure (described in Section |4.1[) en- 
sures that these objects provide only a negligi- 
ble contribution to the count of detected TN SNe. 
Three other objects (f0749t, f0848r, SN2001fn) 
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Table 1 

Redshift Validation Subset : 23 Spectroscopically Observed TN SNe 



Name RA(J2000) DEC(J2000) 7 TN z spec z SOF T a 



SN2001fo 


04 


:37: 


41 


.45 


-01:29:33.1 


1 


.000 





.772 





.716 


+0 
-0 


033 
.071 


SN2001fs 


04 


:39: 


30 


.68 


-01:28:21.9 





.928 





.874 


1 


.169 


+0 
-0 


085 
.085 


SN2001hs 


04 


:39: 


:22 


.39 


-01:32:51.4 





.998 





.833 





.946 


+0 
-0 


137 
110 


SN2001hu 


07 


:50: 


35 


.90 


+09:58:14.2 





.991 





.882 


1 


.048 


+0 
-0 


085 
.175 


SN2001hx 


08 


:49: 


24 


.61 


+44:02:22.4 





.999 





.799 





.899 


+0 
-0 


123 

.059 


SN2001hy 


08 


:49: 


:45 


,85 


+44:15:31.8 





.998 





.812 





.886 


+0 
-0 


149 
.072 


SN2001iv 


07 


:50: 


13 


.53 


+10:17:10.4 


1 


.000 





.397 





.410 


+0 

-(1 


019 
.084 


SN2001iw 


07 


:50: 


39 


.32 


+10:20:19.1 


1 


.000 





.340 





.391 


+0 
-0 


030 
.054 


SN2001ix 


10 


:52: 


48 


.92 


+57:07:29.6 





.921 





.711 





.642 


+0 

-(1 


098 
.059 


OATOAA1 * 

SN2001iy 


10 


52 


24 


28 


+57:16:36.1 


1 


.000 





.568 





.591 


+(1 

-0 


175 
.045 


SN2001jb 


02 


;2G: 


33 


.31 


+00:25:35.0 





.936 





.698 





.811 


+0 
-0 


111 

.084 


SN2001jf 


02 


:28: 


:07 


.13 


+00:26:45.1 





.881 





.815 





.851 


+0 
-0 


123 
.072 


SN2001jh 


02 


:29: 


:00 


.29 


+00:20:44.2 


1 


.000 





.885 





.764 


+0 
-0 


033 
.072 


SN2001jm 


04 


:39: 


13 


.82 


-01:23:18.2 





.999 





.978 





.883 


+0 
-0 


058 
.072 


SN2001jn 


04 


:40: 


12 


.00 


-01:17:45.9 





.897 





.645 





.797 


+0 
-0 


045 
.189 


SN2001jp 


08 


:46: 


31 


.40 


+44:03:56.6 





.640 





.528 





.512 


+0 
-0 


071 
.046 


SN2001kd 


07 


:50: 


31 


.24 


+10:21:07.3 





.934 





.936 





.893 


+0 
-0 


111 

.059 


SN2002P 


02 


:29: 


:05 


.71 


+00:47:20.1 





.772 





.719 





.573 


+0 
-0 


137 

.098 


SN2002W 


08 


:47: 


:54 


.42 


+44:13:42.9 





.822 


1 


.031 


1 


.067 


+0 
-0 


111 

149 


SN2002X 


08 


:48: 


30 


,54 


+44:15:35.3 





.939 





.859 


1 


.012 


+0 
-0 


071 

163 


SN2002aa 


07 


:48: 


:45 


.28 


+ 10:18:00.8 





.926 





.946 





.826 


+0 
-0 


123 

.045 


SN2002ab 


07 


:48: 


55 


,70 


+ 10:06:06.3 





.993 





.423 





.471 


+0 
-0 


033 
111 


SN2002ad 


10 


:50: 


12 


,19 


+57:31:11.6 





.996 





.514 





.484 


+0 
-0 


163 

.045 



a The z estimate determined by SOFT without a spectroscopic redshift prior. 
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Table 2 

SOFT results for 131 TN SN Candidates 



Name 


RA(J2000) 


DEC(J2000) 


7tn 




zsoft 3. 


z .b 


Zprior Source 


olNzUUln 


UZ.ZO.OO.O ( 


-f-UU.ZD.lO.I 


1 


00 


u 


33 


+0.01 
-0.01 


U.ozU ± U.UUl 


T A TTf"' 77 A K 
LAV 1 1 40 


SN2001£k 


02:28:04.70 


+00:46:17.7 


1 


00 





81 


+0.01 
-0.01 


0.804 ± 0.001 


DEEP2 


SN2001fn 


02:29:00.49 


+00:42:21.6 





69 





76 


+0.02 
-0.04 


0.186 ±0.01 


IAUC 7745 


SN2001fo 


04:37:41.45 


-01:29:33.1 


1 


00 





77 


+0.01 
-0.01 


0.772 ± 0.001 


B04 


SN2001fq 


04:38:04.51 


-01:23:49.0 


1 


00 





94 


+0.18 
-0.06 


0.01+0.01 


IAUC 7745 


f0230a 


02:26:09.11 


+00:34:56.9 


1 


00 





77 


+0.01 
-0.01 


0.771 ± 0.001 


DEEP2 


f0230b 


02:27:26.13 


+00:48:27.7 





53 





24 


+0.00 
-0.02 






f0230c 


02:27:54.90 


+00:31:48.4 


1 


00 


1 


00 


+0.06 
-0.10 






f0230e 


02:28:21.92 


+00:24:26.8 





60 


1 


02 


+0.01 
-0.01 


1.014 ± 0.001 


DEEP2 


f0230f 


02:29:35.86 


+00:38:45.2 





05 





86 


+0.08 
-0.06 







* The complete table is available in the online version. 

a Redshift determined with SOFT using the best available z prior 

b Independent SN or host galaxy redshift measurement, where available. 

c Source of the redshift prior: 
IAUC = Spectroscopic z measurement reported in IAU Circulars; 

DEEP2 = Spectrocopic z from the DEEP2 Redshift Survey, Data Release 2 

( |http : //deep . berkeley . edu/DR2[ ); 

SDSS = Photo-z from the SDSS Data Release 6 photoz2 catalog (jOvaizu et al.ll2008l ). 
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have SOFT redshift estimates that are signifi- 
cantly higher than their priors. The priors for the 
first two are again from SDSS galaxies, and in this 
case both galaxies have > 5" separation from the 
position of the SN. We assume that these are fore- 
ground galaxies, and the actual hosts have a low 
surface brightness and arc undetected by SDSS. 
The final object, SN2001fn, is well aligned with a 
galaxy in the SDSS that has a photometric redshift 
of 0.7±0.3, although the reported spectroscopic 
redshift of the host galaxy is 0.186. Even with a 
strong spectroscopic prior at z=0.186, SOFT finds 
a good classification for this object as a TN SN at 
a redshift of 0.761[J;q4 - in good agreement with 
the SDSS photo-z. Without conclusive evidence 
in either direction, we adopt the SOFT redshift 
estimate, noting that a change in the redshift of 
this single object does not significantly affect our 
final rate measurements. 

4. The Revised IfA Deep TN SN Rate 

The fundamental measurement needed to com- 
pute the TNSN rate as a function of redshift is 
simply a count of the number of detected SNe 
as a function of redshift: N^ ct (z). This count 
is computed in Section 14.11 for 9 redshift inter- 
vals reaching out to z=l.l. To convert this to a 
volumetric rate, the detection count at each red- 
shift is divided by the control count, N ct -ri(z) = 
SNuVol x V(z) x i c triO). Here V(z) is the vol- 
ume of the survey in the redshift range from z 
to z + dz and t c tri(z) is the effective survey pe- 
riod or control time: the total time span in which 
a TN SN exploding at redshift z would be de- 
tectable by this survey. Thus, the control count 
gives the number of SN detections that would 
be expected from this survey if the SN rate per 
unit volume were constant at all rcdshifts, fixed 
at 1 SNuVol=lCT 4 SNIa Mpc~ 3 yr" 1 h% . Simu- 
lations to compute N c tvi(z) are described in Sec- 
tion [O 



4.1. The Detected Count 

In the absence of a strong spectroscopic red- 
shift prior, the redshift distributions produced by 
SOFT are typically broad and asymmetric, some- 
times with multiple peaks - making it impossi- 
ble to pin them down into any single redshift 
bin for counting. To incorporate these complex 
uncertainties into our TNSN detection counts, a 



Markov Chain Monte Carlo (MCMC) sampling 
test was executed for each SN, as illustrated in 
Figure H The MCMC provides 1000 redshift re- 
alizations for each SN, distributed so as to recon- 
struct the measured SOFT redshift distribution. 

After repeating this for all candidates, we have 
1000 sets of 131 redshift assignments, with each set 
representing a possible redshift distribution for the 
entire IFA Deep SN sample. To count the num- 
ber of detected SNe, each set of 131 redshifts is 
binned into intervals of width Az=0.1 from 0.0 
to 1.4. To incorporate the uncertainty of clas- 
sification into our counts, we allow each SN to 
contribute a fractional count to the redshift bin 
that it falls in for a given simulation. This frac- 
tional count is drawn from a Gaussian distribu- 
tion centered on the composite TNSN classifica- 
tion grade, 7tn, with a width of a = 0.05. This 
5% scatter reflects the uncertainty of the classifi- 
cation grades ded uced from testing SOFT against 
validation data in iRodnev fc Tonrvl (|2009l) . With 
each SN contributing a count of less than unity, 
the average of the total count across all redshifts 
is S7tn = 84.5, rather than 131. 

After binning up all 1000 Monte Carlo realiza- 
tions, we have 1000 measurements for the count 
in each bin. Taking the mean of this 1000-clcment 
vector gives us our central count estimate for that 
redshift interval, and the standard deviation pro- 
vides an estimate of the statistical uncertainty 
that accounts for error in both the classification 
and the redshift estimation. The counts of de- 
tected SNe in each redshift bin are given in Ta- 
ble [5^,, along with the statistical and systematic 
uncertainties (described in Section [SJ. 

4.2. The Control Count 

The second component of the TN SN rate cal- 
culation, iVctri, is derived from a Monte Carlo sim- 
ulation of ^5000 SNe that replicates the IfA Deep 
survey observations. The simulation is defined 
over a grid of redshift values from z = 0.01 to 1.6 
in steps of 5z = 0.02. At each point on the grid we 
compute -/V oxp (;z), the total number of SNe that ex- 
plode each year within the IfA Deep survey area in 
a shell of thickness 5z. This number is calculated 
as N cxp {z) = 1 SNuVol x V(z) x 1 yr. Here V{z) 
is the volume of a shell at redshift z, using a total 
survey area of 2.29 square degrees and assuming a 



6 



Table 3 

New Determination of IfA Deep TN SN Rates 



z a 


iVdot 


Detected Count 
5N Z SN Poiss 


b 

(SN syst ) 


Control Count c 
iVotrl (5N syst ) 


SNR 


SN Rate 
SSNR stat 


d 

(5SNR S!/St ) 


0.15 


1.95 


±0.12 


±1.40 


(+0.04\ 
\-0.32) 


6.1 


H: 6 3 ) 


0.32 


±0.23 


/+0.07\ 
l-0.06j 


0.35 


4.01 


±0.91 


±2.00 


(+0.29\ 

\-o.oo) 


11.7 


(tli) 


0.34 


±0.19 


/+o.07\ 

1-0.03,/ 


0.45 


5.11 


±1.14 


±2.26 


(+l.74\ 
l-O.W 


16.6 


( + -f 2 o) 


0.31 


±0.15 


(+0.12\ 
\-0.04) 


0.55 


6.49 


±1.42 


±2.55 


(+0.17\ 
\-1.30) 


20.4 


(tli) 


0.32 


±0.14 


(+0.07\ 
\-0.07) 


0.65 


10.09 


±1.75 


±3.18 


1,-1.09 ) 


20.8 


m) 


0.49 


±0.17 


(+0.U\ 
1-0.08J 


0.75 


14.29 


±2.19 


±3.78 


/+0.60\ 
1-2. 03} 


20.9 




0.68 


±0.21 


/+0.23N 
1-0.14J 


0.85 


15.43 


±2.09 


±3.93 


H-1.78N 

l-o.ooj 


19.9 


(±#) 


0.78 


±0.22 


(+a.3i\ 

\-0.16) 


0.95 


13.21 


±2.31 


±3.63 


/+0.32\ 
1-1.43 ) 


17.3 


(± 5 7 :D 


0.76 


±0.25 


(+0.32\ 
l-0.26j 


1.05 


11.01 


±2.08 


±3.32 


/+1.28\ 
\-l.S3l 


13.9 




0.79 


±0.28 


/+0.36N 
\-0Al) 



a Central redshift for bins of width 0.1 (except for the first bin, centered on z = 0.15, which has width 0.3) 

k-?Vdet gives the number of detected TNSNe in each bin, and SN Z is the uncertainty due to the redshift 
distributions, described in Section 14.11 The Poisson uncertainty SNp [ BS is simply V^doti an( l the systematic 
uncertainty SN ByB t is taken from the Table [4^. 

c A r pre( j gives the control count of TN SNe as derived in Scction l4.2l The systematic uncertainty 8N ByB t is taken 
from Table gp. 

^The derived TNSN rate measurements in units of 10 -4 SNIa Mpc -3 yr _1 h~ , computed as SNR= 
A'det/A'ctri, with accompanying statistical (<5SNR s t a t) and systematic (5SNR sys t) uncertainties. 
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Fig. 2. — To tally the contribution from each individual SN candidate over a series of redshift bins, we utilize 
a Markov Chain Monte Carlo series, drawing 1000 random redshifts that reproduce the redshift distribution 
measured by SOFT, (a) The composite fuzzy membership function (equivalent of a posterior probability 
distribution) for SNf0230a in the (z,/j. e ) plane, (b) After marginalizing over fi e , the distribution in redshift 
(solid line) is used as the target distribution for our MCMC series. The blue histogram shows the results of 
1000 MCMC draws, effectively reproducing the highly irregular probability distribution. 



ACDM cosmologyH 

Having computed the total number of explo- 
sions for each point on the redshift grid, we now 
want to simulate how these SNe would be dis- 
tributed in magnitude space. We assume a lu- 
minosity function for the TN SN population com- 
posed of three Gaussian distributions representing 
the overluminous objects (Ia+ or SN 1991T-like), 
which make up 20% of the population, normal la's 
at 64%, and underluminous TN SNe (la- or SN 
1991bg-like) comprising 16% (|Li et al.ll200lh . The 
Gaussians have a =0.45, 0.30, and 0.50 magni- 
tudes for the Ia+, la, and la— components, re- 
spectively. The same luminosity function is used 
for all redshift bins. 

Each of the ^5000 simulated SNe now has an 
assumed redshift z and a peak magnitude rn p k 
drawn from the luminosity functions. These pa- 
rameters are then applied to construct a model 
light curve using the apparatus of the SOFT 
template library (jRodnev fc Tonrvl 120091) . The 



simulated light curve includes cosmological dim- 
ming appropriate for the redshift by assum- 
ing a ACDM cosmology, as well as extinction 
due to host galaxy dust that is drawn from a 



gala ctic line of sight" di stribution (|Hatano et al 



1998) as implemented in Tonrv et al.l ( 2003 ) and 



Wood-Vasev et al.l (|2007l ) . Briefly, we draw a host 
extinction value for each simulated SN from a 
probability distribution constructed by assuming 
that 25% of hosts are bulge systems and 75% 
are disk systems. The exti nction distribution 
curves of lHatano et al.l (1998) are approximated 
by f(A B ) cx 0.02S(A B ) + lO" 1 - 25 "^ for bulge 
systems and f{A B ) oc 0.02S(A B ) + e -3-77-A 2 B for 
disk systems. 

The control time i c t r ;( z i m pk) for this model 
light curve is calculated by counting the days on 
which this SN could explode and subsequently be 
detected by the IfA Deep survey. To determine if 
a model SN would be detected, we count up the 
number of IfA Deep observation epochs for which 
the synthetic light curve is brighter than the de- 



3 A flat, A-dominated univers e with Ho=71, f!jyj=0.26, 
Q A =0.74 l lKomatsu et alj|2010ft . 
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tection threshold^ If the tally of detection epochs 
is > 5 then we consider this SN to be detected. 

Multiplying together the number of SN explo- 
sions N exp (z,m p k) per year and the control time 
tctri(z,rn p k) gives us a final estimate for the con- 
trol count of SN detections N ctr i{z,m p k)- Inte- 
grating over all possible peak magnitudes gives us 
the final control count of TN SN detections in each 
redshift bin, presented in Table [5Jd. The system- 
atic uncertainty estimates included in Table [3] arc 
derived in the following section. 

5. Systematic Uncertainty Estimates 

Both of the principle components of our rate 
calculation - the detection count and the control 
count - can potentially introduce systematic er- 
rors into the SN rate values. A suite of systematic 
uncertainty tests were performed in which a pa- 
rameter of interest was adjusted, the classification 
and redshift estimation were repeated for all can- 
didates, and new rate estimates were measured. 
The difference between the revised rate and our 
baseline rate measurement was taken as the sys- 
tematic error estimate, separately for each redshift 
bin. All tests are described below, and the result- 
ing uncertainty estimates are reported in Table |U 

5.1. Systematic Errors in N^ et 

Peculiar SNe: If a normal TNSN gets 
matched by a highly peculiar SN template we 
could derive an erroneous redshift. We don't ex- 
pect a significant number of peculiar la's in a 
survey of this size. This systematic test removes 
the very unusual SN 2000cx and SN2002cx from 
the SOFT template library. 

Host Galaxy Dust: First, we change the 
ratio of total to selective extinction to Ry=2.1, 
accommodating the possibility that dust around 
extragalactic SNe is very different fro m the Milky 
Way extinction law of Ry=3.1 (lAstier et a\ 



2006; Kowalski et al 



20081 iKessler et al.l l2009t 



Poznanski et al.ll2009l ). Next, we change the prior 
on the extinction parameter Ay from the "Galac- 
tic Line of Sight" prior to an exponential distri- 
bution p(Ay) oc exp(— Ay/5) tha t approximates 



the extreme inclination models of Riello &: Patat 



(120051) . following iNeill et al.1 (120061) . 



Redshift Prior: For objects without a spec- 
troscopic redshift measurement or a host galaxy 
photo- z, we have relied on a redshift prior that 
scales with the cosmological volume clement, go- 
ing as z 2 . This can introduce a systematic bias to- 
wards higher z estimates. An alternative and dif- 
ferently biased redshift prior can be derived from 
the control count, as given in Table [3] The control 
count vs. z for this survey is well approximated 
by a Gaussian with /i = 0.7 and a = 0.38, so we 
estimate the systematic uncertainty by using this 
function for the redshift prior. 

"Fuz zy AND" Rejection T hreshold As de- 
tailed in lRodnev fc Tonrvi (|2010l) . the SOFT red- 



4 As was done for the actual survey, the I band is the only 
search filter, and the detection limit is set at the single- 
observation 5-cr point source limit of 1=25.2 mags. 



shift estimate for any candidate SN is derived by 
combining the fuzzy membership functions from 
multiple templates using the Fuzzy AND opera- 
tion. To prevent the combination from being dom- 
inated by extreme outliers, we have used a rejec- 
tion threshold of 5% on the integrated membership 
grade jtn- To test whether this value is introduc- 
ing systematic effects, we reprocessed the SN light 
curves using threshold values of 3% and 10%. 

5.2. Systematic Errors in A^ ctr i 

TN SN Luminosity Function: There is sub- 
stantial uncertainty in the balance of the SN lumi- 
nosity function between the overluminous Ia+, the 
normal la, a nd the underlum inous la— sub-classes, 
provided by iLi et "ah ( 200ll ). We first modify the 
baseline (20%, 64%, 16%) to (13%, 77%, 10%) and 
then to (27%, 51%, 22%). 

Detection Thresholds: The actual process 
of extracting 131 likely TNSN candidates from 
among the many thousands of potential transient 
objects in the I fA Deep surv ey was quite complex 
(for details see Barrisl ( 20041 )). Our proxy for the 
selection process used the simple requirement that 
a SN must be above the 5-c point source threshold 
in 5 epochs in order to be counted. Our system- 
atic test modifies this to 4 and then 6 epochs. The 
reported systematic error estimate is half of the 
difference between the baseline and the modified 
control count. 

Survey Area: The computation of the survey 
area was greatly simplified in our Monte Carlo pro- 
cedure by assuming a 90% filling factor for volume 
calculations. We recomputed the control count us- 
ing alternate filling factors of 86% and 94%. 
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Table 4 

Systematic Uncertainty Estimates 



Redshift 


: 0.25 


0.35 


(a 

0.45 


) Change in Detected Count 
0.55 0.65 0.75 0.85 


0.95 


1.05 


drop peculiar templates 


-0 


00 


+0 


,02 


+0 


,72 


+0 


00 


-0.69 


+0.22 


+0 


.06 


— 1.33 


+0.93 


change Ry by 1 


-0 


.32 


+0 


.14 


+ 1 


,42 


+0 


17 


-0.72 


-0.89 


+0 


.08 


-0.34 


+0.20 


flat Ay prior 


+0 


.04 


+0 


,20 


+0 


,71 


-0 


32 


+0.00 


-0.50 


+0 


,71 


-0.39 


-0.77 


Gaussian z prior 


+0 


.01 


+0 


.10 


-0 


,06 


-0, 


08 


+0.34 


-0.52 


+0 


13 


-0.02 


-0.02 


modified fuzzy AND threshold . . 


+0 


00 


+0 


.11 


-0 


.15 


-1. 


26 


+0.89 
-0.46 


+0.56 
-1.67 


+1 


.63 


+0.32 
-0.04 


+0.85 
-1.66 


TOTAL : 


+0 


.04 


+0 


.29 


+1 


,74 


+0 


17 


+0.95 


+0.60 


+1 


,78 


+0.32 


+1.28 


-0 


.32 


-0 


.00 


-0 


.16 


-1. 


30 


-1.09 


-2.03 


-0 


.00 


-1.43 


-1.83 



Redshift: 



0.25 



0.35 



(b) Change in Control Count 
0.45 0.55 0.65 0.75 0.85 



0.95 



1.05 



Alternate SN luminosity function 



+0.2 
-0.3 



+0.4 
-0.5 



+1.0 +0.4 



-0.5 



-1.1 



+0.5 
-0.8 



+0.0 
-0.7 



+0.0 
-1.2 



+0.0 
-1.1 



+0.3 
-1.5 



Minimum no. detections ± 1 



+0.55 +1.1 
-1.05 -1.8 



+ 1.95 +2.3 
-2.25 -2.6 



+2.45 +3.05 +4.1 
-3.15 -4.45 -5.6 



+5.65 +6.8 
-5.7 -5.15 



Filling factor ± 4% 



+0.00 +0.05 +0.15 +0.0 +0.0 +0.0 +0.0 +0.0 +0.0 
-0.35 -0.55 -0.75 -1.5 -1.45 -1.85 -2.05 -1.65 -1.35 



Exponential dust profile 



-0.65 



-1.1 



-1.7 



-2.9 



-4.2 



-5.1 



-4.8 



-3.8 



-2.5 



TOTAL : 



+0.59 +1.17 +2.20 +2.33 +2.50 +3.05 +4.10 +5.65 +6.81 
-1.32 -2.24 -2.96 -4.32 -5.47 -7.02 -7.75 -7.13 -6.07 
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> Botticella et al. 2008 ■ Kuznetsova et al. 2008 
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O Neffl et al. 2006 it Tonry et al. 2003 

# Pain et al. 2002 A Barris and Tonry 2006 
O Poznanski et al. 2007 V tni s work 

• Strolger et al. 2003 




1.0 

Redshift 



Fig. 3. — Results of the final TN SN rate calculation. The orig inal IFA Deep rates (jBT06l ) are shown 
as upward-pointing green triangles. The new rate measurements described here are plotted as downward- 
pointing blue triangles. Error bars on the new rate points include both statistical and systematic uncertain- 
tics. A compilation of rates from the literature are plotted as gray symbols, with references as given in the 
legend. 



Host Galaxy Dust: As in the detected count 
systematic uncertainty test above, we modify the 
host galaxy dust first to use Ry = 2.1 and then to 
follow an exponential Ay distribution. 

6. Results 

A summary of the principal components of the 
SN Rate calculation is given in Table El The final 
tally of the number of detected TN SNe is given 
in column 2, with statistical and systematic un- 
certainties in columns 3-5. The control count of 
TNSNe is given in column 6, followed by system- 
atic uncertainty estimates. The resulting SN Rate 
numbers and associated uncertainties appear in 
the final three columns. 

These revised rate estimates are signi ficant!; 
different 



BT0 



from the numbers derived in 
These differences arise from improvements in both 
the measurement of N^et and the calculation of 
N c tri- With an improved SN light curve analy- 
sis program in SOFT, we now have more robust 
photometric redshift estimates and classification 
grades for all 131 SNe. SOFT also has new access 



to more detailed prior information with the addi- 
tion of 36 host galaxy redshift measurements from 
spectra and photo-z's. These improvements have 
allowed us to reliably extend the r ate est imate be- 
yond the z=0.75 limit applied by IBT06L reaching 
now out to a mean redshift of 1.05 in the final bin. 

Using an updated Monte Carlo simulation, we 
have generated a revised control count for the 
number of SNe detectable by the IfA Deep survey 
as a function of reds hift. O ur simulations consis- 
tently disagree with lBT06l in that we compute a 
larger control count in every one of the redshift 
bins. However, we have also performed a scries 
of systematic uncertainty tests, and we find that 
when systematic errors are included, our co ntrol 
counts are consistent with those from lBT06l 

In Figure[3]we show the final rate measurements 
from the IfA Deep survey, compared to other rates 
from the literature. The original IfA Deep rates 
are shown in green, and the newly derived rates 
are plotted in blue. A compilation of rates from 
the literature is plotted in the background in gray. 
The new IfA rate measures are clearly in much 
better agreement with the consensus of literature 
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rates than were the IBT06I values. Our new IfA 



Pain, R., et al. 2002, ApJ, 577, 120 



results indicate that the TN SN rate increases with 
redshift out to at least z=1.0, with no evidence of 
reaching a peak in the redshift range sampled. 



— . 2007, MNRAS, 382, 1169 



Poznanski, D., ct al. 2009, ApJ, 694, 1067 
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